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Abstract 

We consider K links operating concurrently in the same spectral band. Each 
transmitter has multiple antennas, while each receiver uses a single antenna. 
This setting corresponds to the multiple-input single-output interference chan- 
nel. We assume perfect channel state information at the single-user decoding 
receivers whereas the transmitters only have estimates of the true channels. The 
channel estimation errors are assumed to be bounded in elliptical regions whose 
geometry is known at the transmitters. Robust beamforming optimizes worst- 
case received power gains, and a Pareto optimal point is a worst-case achievable 
rate tuple from which it is impossible to increase a link's performance without 
degrading the performance of another. We characterize the robust beamform- 
ing vectors necessary to operate at any Pareto optimal point. Moreover, these 
beamforming vectors are parameterized by K{K — 1) real-valued parameters. 
We analyze the system's spectral efficiency at high and low signal-to-noise ratio 
(SNR). Zero forcing transmission achieves full multiplexing gain at high SNR 
only if the estimation errors scale linearly with inverse SNR. If the errors are 
SNR independent, then single-user transmission is optimal at high SNR. At 
low SNR, robust maximum ratio transmission optimizes the minimum energy 
per bit for reliable communication. Numerical simulations illustrate the gained 
theoretical results. 
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1. Introduction 



We consider multiple transmitter-receiver pairs simultaneously operating on 
the same frequency band. The signal from a transmitter is useful information at 
the intended receiver while it is regarded as interference at unintended receivers. 
All transmitters are equipped with multiple antennas while the receivers use only 
a single antenna. This setting corresponds to the multiple-input single-output 
(MISO) interference channel (IFC) Ij. 

In the MISO IFC, the beamforming vectors used at the transmitters impact 
the performance of the systems. A jointly efficient operating point corresponds 
to a Pareto optimal point in which it is not possible to improve the performance 
of one link without degrading the performance of at least another link. Designing 
a Pareto optimal mechanism requires finding the joint beamforming vectors 
used at the transmitters that lead to the Pareto optimal point. In the MISO 
IFC, finding specific Pareto optimal points such as the maximum sum-rate or 
proportional-fair is proven to be strongly NP-hard The importance of 

characterizing the set of beamforming vectors necessary for the links' Pareto 
optimal operation is twofold. First, the set of relevant beamforming vectors 
is reduced to a considerably small subset of all feasible beamforming vectors. 
Second, the characterized set of efficient beamforming vectors is parameterized 
by a number of scalars which can even reduce the complexity of indicating the 
required beamforming vectors. 

Characterizing the beamforming vectors necessary to achieve all Pareto op- 
timal points in the MISO IFC has been carried out in several works for the 
case of perfect channel state information (CSI) at the transmitters. In [ij, the 
efficient beamforming vectors are parameterized by K{K — 1) complex- valued 
parameters, where K is the number of links. For the special two-user case, the 
efficient beamforming vectors are proven to be a linear combination of max- 
imum ratio transmission and zero forcing transmission. The extension to a 
real- valued parametrization for the general if-user case is conducted in 0, @, 0| 
where K{K — 1) real- valued parameters are required to achieve all Pareto opti- 
mal points. Recently in Q, parametrization of the efficient beamforming vector 
is provided in the multi-cell MISO setting with general linear transmit power 
constraints at the transmitters. For the case of MISO IFC and total power 
constraint at the transmitter, the number of required parameters is 2K — 1. In 
[il, : the beamforming vectors necessary and sufficient to achieve all Pareto 
optimal points are characterized and parameterized by a single real-valued pa- 
rameter for the two-user MISO IFC case. If the transmitters only know the 
covariance matrices of the channel vectors, characterization of the Pareto op- 
timal transmission strategies is done for the two- user MISO IFC in [ll|. In 
a multi-cell MISO setting with uncertainty in CSI at the transmitters, robust 
Pareto optimal beamforming is obtained by robust fairness-profile optimization 



12| . All Pareto optimal points in the performance region are achieved re- 



quiring K — 1 real-valued parameters. In addition, a monotonic optimization 
algorithm is applied to achieve specific Pareto optimal points such as maximum 
sum and proportional fair performance points. 
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The presence of CSI at a transmitter is essential in order to increase the 
performance of the multi-antenna system through sophisticated beamforming 
techniques [l3| . In practical communication systems, CSI at the transmitter is 
usually not perfect |14| . The transmitter gains CSI either through reciprocity of 
the uplink and downlink channels or through a feedback link from the receiver 



15| . Both techniques of acquiring CSI at the transmitter entail a delay which 



leads to outdated channel information and hence to a mismatch to the true 
channels. In the feedback model, the receiver quantizes the channel vector 
according to a vector codebook and sends the transmitter the channel index 
which is nearest to the true channel. This mechanism reduces the feedback bits 
required to identify the channel at the transmitter [l^ . The size of the codebook 
used determines the feedback overhead and the accuracy of the selected channel 



vector to the true channel 17 



With imperfect CSI at the transmitter in a multi-user MISO downlink sys- 
tem, robust beamforming is studied in [l^ to minimize the worst-case mean- 
square-error (MSE), and in [l9l.[20| the problem of minimizing the transmission 
power subject to worst-case quality-of-service (QoS) constraints at the receivers 
is analyzed. Furthermore, robust transceiver design in a multi-user multiple- 
input multiple-output (MIMO) system is considered in 2l|,|22|. In the multi-cell 
multiuser MISO setting with imperfect CSI at the transmitters, the problems 
of maximizing the worst-case weighted sum rate and minimizing the weighted 
sum transmission power subject to worst-case QoS constraints are studied in 



23| and [2j|, respectively. 

In this paper, we assume that the transmitters have imperfect CSI. We 
adopt a deterministic uncertainty model in which the channel estimation error 
is bounded in an uncertainty region. The channel uncertainty region is assumed 
to be an ellipsoid [2^ whose geometry is known at the transmitter. Robust 
beamforming takes into account the worst-case channel estimate in the un- 
certainty region. We derive the worst-case power gains at the receivers and 
formulate accordingly the worst-case achievable rates for the links. The robust 
rate region is the set of all worst-case jointly achievable rate tuples. In order to 
characterize the Pareto optimal points of the robust rate region, we show that 
the objective of each transmitter is a tradeoff between maximizing the intended 
worst-case power gain and minimizing the worst-case interference gains. This 
multi-objective problem is cast as a second order cone program (SOCP) which 
can be solved efficiently. Consequently, we characterize the beamforming vec- 
tors that are necessary to achieve all Pareto optimal points in the robust rate 
region. Moreover, these beamforming vectors are parameterized by K{K — 1) 
real- valued parameters taking values between zero and one. Afterwards, we ana- 
lyze the spectral efficiency of the system for asymptotic values of signal-to-noise 
ratio (SNR). At high SNR, achieving full multiplexing gain with zero forcing 
transmission requires the channel estimation error to reduce linearly with SNR. 
If the error is independent of SNR, then single-user transmission is optimal. 
At low SNR, it is shown that joint maximum ratio transmission is optimal to 
minimize the minimum energy per bit for reliable communication. 

The paper is organized as follows. In Section [21 we describe the system 
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model as well as the channel uncertainty model. The worst-case power gains 
at the receivers are derived and the worst-case achievable rates for the links 
are formulated. In Section [3J we characterize the beamforming vectors that 
are necessary to achieve Pareto optimal points in the robust rate region. As 
a special case, the two-user MISO IFC with spherical channel uncertainty is 
analyzed. In Section 21 spectral efficiency in the high and low SNR regime are 
studied. In Section [5l we draw the conclusions. 

Notations 

Column vectors and matrices are given in lowercase and uppercase bold- 
face letters, respectively. ||a|| is the Euclidean norm of a G C''^. |&| is the 
absolute value of 6 € C. {■)'^ and (•)^ denote transpose and Hermitian trans- 
pose, respectively. The orthogonal projector onto the column space of Z is 
Hz := Z{Z^ Z)'^ . The orthogonal projector onto the orthogonal comple- 
ment of the column space of Z is 11^ .— I Tlz, where I is an identity matrix. 
CA/'(0, A) denotes a circularly-symmetric Gaussian complex random vector with 
covariance matrix A. (a)+ denotes max(a, 0). Re{a), Im{a) and /(a) denote 
the real part, imaginary part, and the phase of a complex number a. 



2. System Model 

Consider a ivT-user MISO IFC. Each transmitter has an intended receiver, 
and the signal from a transmitter to an unintended receiver is treated as inter- 
ference. Each transmitter k is equipped with Nk antennas while all receivers use 
single antennas. The quasi-static block flat-fading channel vector from trans- 
mitter k to receiver i is denoted by h^i G C^*". We assume that transmission 
consists of scalar coding followed by beamforming. The beamforming vector 
used by transmitter k is Wk S C^''. Each transmitter k has a total power 
constraint of Pk such that Hii'fcp < Pk- The matched-filtered, symbol-sampled 
complex baseband data received at receiver £ is 

yi = 2^^^^hk(.WkSk+ni, (1) 

where Sk ~ CA/'(0, 1) is the symbol transmitted by transmitter k. The random 
variables Uk ^ CJ\f{0,a^) are additive Gaussian noise. 

We assume that perfect CSI is present at the receivers while at the trans- 
mitters the channels are not perfectly known. A transmitter has estimates of 
the true channels between itself and all receivers, and we assume that the error 
in the channel estimation vectors is bounded in an uncertainty region. The 
ellipsoidal channel uncertainty region is described next. 

2.1. Channel Uncertainty Model 

Let hki be the estimate of the true channel vector hki at transmitter k. The 
uncertainty in the channel estimate can be modeled using a channel estimation 
error vector 27 1 

<5fc£ — hki — hkt, Ski G Eki, (2) 
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where Ske is assumed to be bounded in an ellipsoidal uncertainty region £ki 
defined as 

Ski = {Am5 : \\5\\ < em}. (3) 

In ([2]), Aki S C^k'^N^ determines the shape of the elhpsoid. We assume that 
£ki is full rank for all k,£ — 1, . . . , K, and has a maximum radius of eki, i-e. the 
largest singular value of A^e is constrained to one. Generally, the shape and size 
of the ellipsoid in ([3]) should be specified according to the error probability in 
channel estimation for the system under consideration. For example, if channel 
vector estimation is done at a transmitter using training-sequences from the 
corresponding receiver, the size of the uncertainty region can be chosen as a 



scaled version of the channel estimation mean square error (MSE) |22l. |12| 



The elliptical uncertainty model described in ([3|) is more general and encom 



passes the spherical uncertainty model adopted in [28|. The spherical uncer- 
tainty region is defined as 

Vki = {d:\\d\\<eke}. (4) 

Specifically, £ke in ([3]) is a sphere when Aki = I- 

A transmitter k knows the channel estimates hke and the associated un- 
certainty region £ke- Robust transmission requires the choice of beamforming 
vectors to be robust to channel estimation errors. Thus, a transmitter has to 
consider worst-case achievable rate at its intended receiver for secure communi- 
cation. 

2.2. Worst- Case Achievable Rates 

The worst-case achievable rate for link £ with single-user decoding is 



Rdw,,...,WK) =log,{l+ , ^."''^'i ), (5) 



where x'^p{wk) is the worst-case signal power from transmitter k to receiver £. 
The worst-case intended power gain x'jp{wi) is the least power gain achievable 
within the uncertainty set £u. That is, 

xii{wi) = min \hf^we\ with hu = hu + du (6a) 



min \huWe + 6uWi\ (6b) 



IhfiWi + d" AuWi\. (6c) 



^In [25j . the ellipsoidal uncertainty region is defined as £ki = {A^eS + Cfcf : < 
where cj.^ defines the center of the ellipsoid S^i- Here, we assume c^i = without loss of 
generality. 



5 



The error vector which minimizes the intended power gain in ^ can be calcu- 
lated following similar steps as in [28']. Using the triangle inequality [2^ 3.2.5], 
we have 



\huWe + S Auwe\>< H ^ .h , • ' 

y\hiiWi\ — \o AggWi] otherwise 

for all S such that \\5\\ < e«. In (O, the lower bound of zero for \hggWi\ < 

{S^Aj^fWil is achievable by a negative linear scaling of (5. Otherwise, if j/i^^iOfl > 
\5^ Af^Wil, then the error vector which achieves the lower bound in (O is 

S = -eee e^Ahf,^.)_ (g) 

Substituting ([5]) in the RHS of ([7]), the square root of the worst-case intended 
power gain in ([6]) reduces to 

xu{w() = (^\h"w£\ - \\AuWe\\eu^ ^. (9) 

The square root of the worst-case interference power gain from transmitter k to 
receiver £,k ^ is 

iki + Ski (10a) 
(10b) 
(10c) 

Using the triangle inequality, we have 

\h",Wk + S"A^,Wk\ < \h",Wk\ + \S"A^,Wkl (11) 
and the upper bound in ()lip is achieved by 

5^ eke ^"r' e^^'-^^'^-l (12) 
Hence, the worst-case interference gain is 

Xke{wk) = \hkiWk\ + \\A"wk\\eke. (13) 

Note that the worst-case measures in © and are the main terms in 
the worst-case rate expression in ([S]). Also, the channel estimation errors of 
all channel vectors are independent such that the worst-case interference and 
intended power gains may occur simultaneously. 



= max 




with hke ~ 






~ max 


\hkiWk H 




Skt££kt 


— max 


\hf^Wk - 


h6"AZwk\. 


\\S\\<eke 
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3. Pareto Optimal Beamforming 



The if-dimensional robust rate region consists of all jointly achievable worst- 
case rate tuples defined as 

n {{Ri{w,,...,Wk),...,Rk{'Wi,...,Wk)) : \\wk\\^ < Pk,k = 1, . . . , K}, 

(14) 

where Rkiwu- ■ ■ , wk) is defined in ([5]). The region TZ in ((Tl)) is not necessarily 
a convex se1o [30|]. This is mainly due to the interference coupling present 
between the links. The outer boundary of the rate region comprises efficient 
operating points. These points are called Pareto optimal. 

Definition 1. A rate tuple (Ri, Rk) £ TZ is Pareto optimal if there is no 
other tuple {R'l, R' k) G such that (i?'i, R' k) ^ (^i, Rk), where the 
inequality is component-wise and strict for at least one component. The set of 
all Pareto optimal operating points constitutes the Pareto boundary ofTZ. 

In order to achieve a Pareto optimal point, each beamforming vector must 
include a tradeoff between the maximization of the intended power gain and the 
minimization of the interference power gains. Next, we will consider the special 
two-user MISO IFC case and compare Pareto efficient beamforming with perfect 
and imperfect CSI at the transmitters. Also as a special case, the channel 
uncertainty regions are considered to be spherical. Later, the general K-usev 
MISO IFC with elliptical channel uncertainty region will be restored. 

3.1. Two-User Case with Spherical Channel Uncertainty 

In the two-user case, the rate region in (1141) is two-dimensional rewritten as 

7^ := {{Ri{wi,W2),R2{wi,W2)) : <Pk,k^ 1,2}, (15) 

where the achievable rate for link £ is 

R,{w,,W2) ^log2(l+ ^^^"""f , A k^L (16) 

If the transmitter knows the channel vectors perfectly, the beamforming vectors 
that are necessary to achieve all Pareto optimal points of the two-user rate 
region in are Q 

Wk{ak) = T^fy^ +^/^^ |in^"!!''A (17) 

where au € [0,<=^^], with = \\nh,,hkk\\' l\\hkk\f . In the beam- 

forming vectors that achieve Pareto optimal points are a combination of two 



''Note that any point on the convex hull of TZ. can be reached by time-sharing. In this work, 
we do not consider time-sharing techniques to convexify the rate region TZ. 
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orthogonal unit norm vectors. Zero forcing (ZF) to the unintended receiver 
corresponds to ak = 0, and maximum ratio transmission (MRT) to ak = Oi"^'^ . 

If the transmitters have imperfect CSI and the channel vector uncertainty 
region is spherical as defined in Q , then the worst-case achievable rate of link 

e is 



Rl{wi,W2) = log2 



ififfwA - ||u;£| 



y (J^ + (\h"/,Wk\ + \\wk\\eke 



2 



k ^ e, (18) 



where the intended and interference power gains are from (|9]) and (fT3|) . respec- 
tively, with Aki = I- According to the rate expression in ([T^ . worst-case signal 
power and interference powers include additive terms influenced only by the 
norm of the beamforming vectors. It is thus expected that robust Pareto opti- 
mal beamforming includes additionally varying transmission power. The Pareto 
boundary of the rate region TZ in (jisp is achieved by the beamforming vectors 



m 



, / ^ nr hkk , ^fefe \ 

Wk{^k,l3k) = V^k[^/Fk ,,J''"r II +^^^ iirT±'r n l^ (19) 

^h.hkkW P 



where ^ e [0, 1], and f3k e [0,/3f with /3f ™- = hkkf /llhkkf ■ The 



parametrization in (|T9)) requires two parameters per transmitter. One parameter 
is to alter the direction of the beamforming vector and one parameter scales the 
transmission power. Interestingly, the structure of the efficient beamforming 
vectors with spherical channel uncertainty in (jl9p is similar to the case of perfect 
CSI at the transmitters in (IT71) . In the channel estimates replace the true 
channel vectors in (|17p. The beamforming vectors corresponding to MRT and 
ZF represent extreme strategies which have the objective of either maximizing 
the power at the intended receiver or minimizing the interference power gain. 
Robust MRT is calculated as 

«;r""= argmax \h"kWk\ - Ww^We^k ^ ^/P'k-^. (20) 

\\'Wk\\<'/K WrikkW 

Thus, to maximize the power gain at the intended receiver in the worst-case 
of spherical uncertainty, the transmitter chooses full power transmission in the 
direction of the estimated channel. In [s^], it is shown that robust MRT cor- 
responds to the channel vector estimate whenever the uncertainty region is 
symmetric. Robust ZF, however is achieved by allocating zero power. This is 
observed in the denominator of (1181) where the interference gain can only be 
zero for \\wk\\ = 0. 

Next, we address the general i^T-user case with general ellipsoidal uncertainty 
regions defined in 1^. 
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3.2. K-User Case with Ellipsoidal Channel Uncertainty 

As in the two-user case, Pareto optimal beamforming in the i^-user case 
requires a tradeoff between maximizing intended power gain and minimizing in- 
terference gains. This is concluded by observing that the worst-case achievable 
rate in ([SJ is monotonically increasing with the direct power gain x'jf^{wi) for 
fixed interference power x'^g{wk), k £. In addition, the worst-case achievable 
rate is monotonically decreasing with the interference power gain for fixed in- 
tended power gain. In other words, increasing the power gain at the intended 
receiver increases its achievable rate and reducing the interference at unintended 
receivers increases the rates of the other links. For a transmitter k, the tradeoff^ 
between the two objectives of increasing x1f.{wk) and reducing x1f{wk), k £, 
can be cast as a multi-objective optimization problem 33, 34] 



Wk = argmax Okiwk) s.t. ||iOfe|p < Pk, 
where the multi-objective function g^. : 0^*= M^^ is defined as 



gke{wk) 



-xlg{wk), k^£; 
'4fc("^fc), k = £. 



(21) 



(22) 



There are several methods to solve the problem in 
weighted sum method 



Tj). One technique is the 



Wkhk) ^ ^^S^^^ Ikkxlk(wk) ~^jkexli{wk) s.t. \\wk\\'^ < Pk, (23) 



where ^ke are nonnegative weights such that X^fci ^ki = 1. In the ii'-user MISO 
IFC with perfect CSI at the transmitters, the beamforming vectors necessary to 
achieve all Pareto optimal points are found in by the optimization problem 
For the perfect CSI case, the objective in (|23p is a sum of Hermitian 



forms and hence the optimization is an eigenvalue problem. 

Another method for solving (|2ip is by maximizing a single objective in 
Qkiwk) while setting goals for the other objectives [s^. This method is called 
the boundary intersection approach [34] and casts the problem in (PT|) as 



Wk{\k) = argmax xlk{wk) 



s.t. 



^kl 



< A^Ffc^, for all £^k, 



kfelr < Pk 



(24a) 
(24b) 
(24c) 



with Afe£ G [0, 1] and P^ is a fixed value corresponding to an upper limit on the 
interference level generated by transmitter k. The formulation in (j24p is common 
for optimization problems in cognitive radio networks where a secondary user 
maximizes his intended power gain or achievable rate subject to interference 

Also, in [6] the 



37 



temperature constraints at the primary receivers '36, 
necessary beamforming vectors to achieve all Pareto optimal points in the MISO 
IFC with perfect CSI are found by maximizing the achievable rate while setting 
interference constraints at the unintended receivers. 
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In order to solve the multi-objective optimization problem in (j2ip with chan- 
nel uncertainty at the transmitters, we use the method in (j24l) and show that it 
leads to a convex optimization problem [39| . By taking the square root of the 
objective function and the constraints and substituting the power gains from 
dH) and (HH) in (IMl) we get 

Wk{\k) = argmax |/iffeiOfc| - WAfkWkWekk (25a) 
s.t. \hkiWk\ + \\AkiWk\\eki < V^M^M, (25b) 

e^i,...,K, i^k, 

\\wk\\ < VK- (25c) 
Since the objective function does not depend on the phase of the beamforming 



vector Wk ^8|, the problem in ((25|) can be equivalently written as 



Wk{\k) = argmax Yle{h"kWk} - WAkkWkhkk (26a) 



s.t. \fh:pWk\ + \\AkiWk\\ekt < V^kt^ki, (26b) 

e^i,...,K, i^k, 

lm{hfkWk} = 0, (26c) 
\\wk\\ < /a, (26d) 

where Xke S [0, 1]. The problem in (f26| is a second order cone program (SOCP) 
(ioj l which can be efficiently solved by interior point methods [39| (An optimiza- 
tion package that solves SOCPs is SeDuMi [41.].) 

In the optimization problem (|26p. the parameter Xke determines the inter- 
ference level that transmitter k is allowed to generate at an unintended receiver 
£. For Xki = 1, the highest interference level is allowed. For the constraints Tke 
with Xk£ = 1 to be tight, Tke should be the interference level at receiver £ when 
transmitter k performs MRT, i.e., 

V.2 R-MRT \ / I I^"^ „..R-MRT I , 1 1 yl -ff „ ,.R-MRT I 



"n ™^ R-MRT\ / II, „ R-MRT I , M y| -H „ R-MRTII \ /0'7^ 

Tkt^XkeiWk }^[\hkeWk \ + \\AkeWk \\eke), (27) 
with MRT beamforming obtained from 

= argmax Re{hkkWk} - WAj^k^kHk (28a) 

s.t. lm{hfkWk} 0, (28b) 

\\wk\\ < /a. (28c) 

Another method for calculating robust MRT is provided in fsil] which uses a 
one dimensional line search. 



Theorem 1. All Pareto optimal points of the robust rate region TZ in (|14p can 
be achieved by the beamforming vectors Wk{X.k) from (j26p with Xk E [0,1]^^^ 
fork^l,...,K. 
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Proof. The proof is by contradiction. Assume a beamforming vector Wk 
which is not a solution of the problem in (j26|) achieves a Pareto optimal point 
(i?i, . . . , Rk)- Then, it is possible to find a beamforming vector loj, from (|26p in 
which the intended power gain increases at receiver k, i.e. x1i^{w'f.) > x\f,{wk)^ 
without affecting the interference gain to all other receivers, i.e. x'^f^{w'^) = 
x1g{wh) for all £ ^ k. In this case the achievable rate of link k is increased 
without affecting the rates of the other receivers. According to Definition [1] 
. . . , Rk) would not be Pareto optimal which is a contradiction to the orig- 
inal assumption. Hence, any beamforming vector not in the solution set of (|26p 
does not achieve a Pareto optimal point in the robust rate region TZ. □ 

The Pareto boundary of the robust rate region can be also obtained by 



the framework provided in [12|. The approach in [12| uses a robust fairness- 
profile optimization to calculate a Pareto optimal point which entails solving 
a set of convex feasibility problems. Accordingly, only points on the Pareto 
boundary are delivered requiring K— 1 real- valued parameters. In comparison to 
the parametrization in [12i | , our result in Theorem [T] specifies the beamforming 
vectors for each transmitter that are necessary to achieve all Pareto optimal 
points in the robust rate region. The efficient beamforming vectors for each 
transmitter are parameterized hy K — 1 real-valued parameters each between 
zero and one. The parametrization of the efficient beamforming vectors can 
be utilized for designing efficient low complexity distributed resource allocation 
schemes which require low signaling overhead between the transmitters as is 
done in [13, ill. 



[Figure 1 about here] 

In Figure [1] a three- user robust rate region is plotted. The total power con- 
straints at the transmitters are set to one and the SNR is 1/cr^ = dB. The 
number of antennas at each transmitter is three. The channel vector estimates 
hke are independent and identically distributed as hki ^ CAf{0,I). Each el- 
lipsoidal uncertainty region is generated as follows: Nk vectors are generated 
according to complex normal distribution with zero mean and covariance matrix 
I. Aki is constructed by concatenating the generated vectors and normalizing 
the constructed matrix such that the largest singular value of Ake is one. The 
errors are equally chosen as ete = 0.5 for k,£ = 1, 2, 3. For each transmitter k, 
we generate the efficient beamforming vectors from ([25)) with the parameters 
Afcf uniformly sampled in a 0.05 step-length between zero and one. A set of rate 
tuples is calculated from the generated beamforming vectors of the transmitters. 
Within this set, a subset corresponds to Pareto optimal points of the robust rate 
region. In Figure (H only the Pareto optimal points are plotted by utilizing the 
MATLAB code in [M]- 

4. Spectral Efficiency at High and Low SNR 

In this section, optimal beamforming is studied for asymptotic values of 
SNR. Moreover, we analyze the effects of imperfect CSI at the transmitters on 
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the performance of the system in comparison to the case of perfect CSI. 



4-1- Efficiency at High SNR 

The quantitative performance is analyzed using the high-SNR offset concept 
in 45|, Section II] . We define the SNR as p = 1 /cr^ and the maximum sum rate 



as a function of SNR as The high-SNR slope is 



Sr, 



p'-^'i log2(p) ' 



lim 



(29) 



which corresponds to the multiplexing gain, i.e. the slope of the maximum sum 
rate curve at high SNR. The maximum sum rate is 



K 



i?=™(p) = max log2 1 + T— 4^ 



(30a) 



K 



= max 



5: log 



p(\h"w(\ - \\Aff>we 



1 + P E (\hfiWk\ + \\A^i,Wk\\ekt 

k^i ^ 



(30b) 



where xu{wi) and Xke{wk) are from ([9]) and p3|) respectively. From ([29]). the 
high-SNR slope is 



max 

wi ^...,wk 



Soo = lim 

p— >oo 




p{\hu-Wi\-\\A"f,wi\\tit^ ^ 
1+P E (\h^t'Wk\ + \\A^t'^k\\^kX 



log2(p) 



(31) 



which is maximized when the interference is nulled at all receivers. Since Am 
is full rank, the worst-case interference gains cannot be nulled unless the trans- 
mitters switch their transmission off. In case the error does not depend on 
the SNR, the interference gains as well as the intended power gain in the achiev- 
able rate of a link I scale linearly with p. Hence, the high-SNR slope in (I5T|) is 
zero if more than one link operate simultaneously. Therefore, in the high-SNR 
regime single-user transmission is optimal achieving the largest high-SNR slope 
of Sao — 1- The maximum sum rate is then 



i?™'(p)=log2 1+P 



max 

= 1 K 



I aH „ .MRT II 

\AuWg ||e« 



(32) 



where only one user operates using MRT and full power transmission. Note that 
the condition that determines the dominant user does not only depend on the 
channel gains but also on the amount of uncertainty present at the transmitter. 

In 17[ it is shown that the channel estimation error has to scale linearly 
with the inverse SNR in order to achieve maximum multiplexing gain with ZF 
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in a MISO broadcast channel. That is, e^eip) °^ p should hold for all k and £. 
We assume that the estimation error has the following dependence on SNR: 

€k£{p) ^ (33) 
where a^e is a constant. ZF beamforming according to the channel estimates is 



where 



hki, . • . , hkk~i, h 



kk+l, ■ 



kK 



(34) 



(35) 



Note that if the number of antennas at a transmitter k is strictly less than the 
number of receivers K , then = 0. Assuming Nk > iiT for all fc = 1, . . . , if, 
the maximum sum rate with ZF is 



R—{p) 



K 

E 

1=1 



K 



log2 




(36a) 



e=i 

K 

£=1 



-H 



^\huwr\-^\\Af,wr\\^^ 



ZF II akt 
k II 



1 



AuwfWau 



The high-SNR slope from (HH) with TVfc > if for aU fc = 1, . . 

X ( {^\ha'^T\-\\Afc-wY\\au 



(36b) 



(36c) 



, K is 



lim ■ 



= if. 



log2(p) ^^^^ 
Hence, the maximum multiplexing gain of K is achieved as with perfect CSI. 
Note, that we assumed that Nk > K for all k. If there exists a transmitter k 
such that Nk < K, then transmitter k cannot perform ZF to K receivers simul- 
taneously. In this case, the highest multiplexing gain would be the maximum 
number of transmitters m* that have more antennas than m* . The optimization 
problem to calculate the multiplexing gain is 

(38a) 
(38b) 



max m 



s.t. 



K 

k=l 



Vk 



1, Nk > m; 
0, Nk < m. 
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Finding m* has to be performed iteratively. In Algorithm [TJ we provide a 
method to calculate m* . First, the vector containing the number of anten- 
nas [Ni, . . . , Nk] is sorted in a weakly decreasing order to [Ni, . . . , Nk] with 
Ni > Ni+i- The multiplexing gain m* is found as the largest k such that 
Nk > k. ZF transmission with Nm'+i antennas reduces the multiplexing gain 
to Nm'+i < Nm* — m* . Consequently, in order to achieve the maximum multi- 
plexing gain, only the links corresponding to the first m* entries in [A^i , . . . , Nk] 
should operate and perform ZF according to the channel estimates. 

[Algorithm 1 about here] 

[Figure 2 about here] 

In Figure [U the maximum sum rate is plotted for two links. Each trans- 
mitter has three antennas and the total power constraints at the transmitters 
are = 1, fc = 1, 2. The channel vector estimates and the uncertainty regions 
are generated as stated in Section 13.21 in the description of Figure [T] The error 
ekj for all k,j = 1,2, is chosen as indicated in the label of Figure [2] where for 
perfect CSI ej^j = 0. For each transmitter fc, we generate the efficient beam- 
forming vectors from (j26p with the parameters Xke uniformly sampled in a 0.001 
step-length between zero and one. The set of rate tuples is calculated from the 
generated beamforming vectors and the maximum sum rate is found by grid 
search. For perfect CSI, full multiplexing gain of two is achieved with joint ZF 
transmission. Also, for the error scaling linearly with the inverse SNR as in (|33)) . 
a multiplexing gain of two is achieved as is obtained in (|37l) . If the error scales 
slower than linearly with inverse SNR, as in ekj = 1/ -y/p, the loss in multiplexing 
gain can be observed. For constant error of ekj — 0.3 the multiplexing gain is 
one where only a single link operates as is derived in (15^ . The transition from 
the operation of the two links to the operation of a single link can be noticed 
at around 30 dB SNR in the maximum sum rate curve for constant error. 

4-2. Efficiency at Low SNR 

In order to study the spectral efficiency in the low SNR regime, two perfor- 
mance measures were introduced in (46| . The first measure is Eb/ Nq^^^ which 
is the minimum energy per bit required for reliable communication. The smaller 
-E-bZ-^Omin fo'" ^ system the more reliable is its operation. The second per- 
formance measure is Sq which is the slope of the spectral efficiency curve at 
^b/No^^^. Since 5*0 describes the growth of the spectral efficiency curve from 
£'f)/-^Omm' larger So is for a system the greater is the gain for increasing en- 
ergy per bit from Eb/No^-^. Here, we use these measures to determine optimal 
beamforming in our setting at low SNR. 

The spectral efficiency Ci{Eb/ Nq^) of link £ is defined as the ratio of trans- 
mission rate Ri to bandwidth B [4^ Section III]: 
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where Eb/No^ is the energy per bit for link i normahzed over the background 
noise. Specificahy, Et/Nog = P£/{NoRe), where Pg is the transmitted power. 
The spectral efficiency function is directly related to the common capacity ex- 
pressioifl dip) through Ce{Eb/NQ^) = dip) for the SNR which solves [46] 



Eb_ 

Nn 



Ciip) 



(40) 



At low SNR, the spectral efficiency function (£'{,/ A^o^) can be expressed as 







(Eb 


_E^ 




~ 3dB 


KNoe 


d.B No^in,e 



dB 



with 



-Bb ^ log, 2 



and So , 



cm 



-cm 



(41) 



(42) 



where and are the first and second derivatives, respectively, with respect 
to SNR. In , Eb I Nq^^^ ^ is the minimum energy per bit required for reliable 
communication for link L It is shown in [47, Proposition 4] that the minimum 
energy per bit in (|42p is also valid in the MISO IFC. Furthermore, it is shown 
that single-user decoding at the receivers in the low SNR regime is optimal to 
minimize Eb/No^-^ ^. The rate of link I from ([5|) is 



=log2 1 + 



(43) 



where xu{wi) and Xkt{wk) are calculated in © and p^ . respectively. For the 
purpose of calculating the minimum energy per bit and the wideband slope in 
(|42|) for a link we have 

Cm^^l(w,l (44) 

and 



Q(0) 



(45) 



From (j44p and (|45p . we can calculate the minimum energy per bit and the 
wideband slope in (H^ as 



Eb_ ^ log, 2 
Nominal xj^iwe) 



and 5*0.1 



2xji,{wi) 



(46) 



The minimum energy per bit from (j46p is minimized by maximizing the intended 
power gain x1^{we), i.e. with MRT beamforming. Thus, joint robust MRT 



^Notice the difference in the notation for the spectral efficiency of linlc I, Ci, and the 
capacity of link t, Ci. 
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beamforming is optimal in the low SNR regime to achieve the minimum energy 
per bit for reliable communication. For perfect CSI at the transmitter 



Eb"''^' log, 2 log, 2 



For imperfect CSI at the transmitters, we have 



Eb'''''' log, 2 



(47) 



(48) 



where lof "^^ is obtained from ([28l)- Comparing (gll) to (|47]), the loss due to 
imperfect CSI can be observed. 

[Figure 3 about here] 

In Figure [31 the Ei,/Nq regions for two links with perfect and imperfect CSI 
are plotted. The E^/Nq region with imperfect CSI is smaller and contained in 
the region with perfect CSI. For perfect and imperfect CSI, joint MRT achieves 
the joint minimum Ei,/Nq. 

The wideband slope from ([^5)) of link £ for perfect and imperfect CSI are 

Sfe'' = — , , (49) 



and 



2(\h"we\ - \\AfgWe\\eu^ 
\hfi,we\ - \\AuWe\\eu) +2J2 (\hkeWk\ + \\A"i,Wk\\e 



respectively. AH jointly achievable wideband slopes for the links constitute a 
slope region [il, |4^ . 

[Figure 4 about here] 

In Figure m the wideband slope regions for two links with perfect and imper- 
fect CSI are plotted. The regions are generated by utilizing the parametrization 
of efficient beamforming vectors in (j26p . The wideband slope region with perfect 
CSI corresponds to the box in which the maximum wideband slope of two for 
both links can be achieved simultaneously (with ZF transmission). The wide- 
band slope region with imperfect CSI is contained in that of perfect CSI, and 
the maximum wideband slope of two is achievable for one link only if the other 
link performs ZF, i.e. switches its transmission off. The wideband slopes with 
joint MRT are marked for perfect and imperfect CSI. These correspond to the 
slope of the spectral efficiency curve at j^^^^^ ^ as is given in (HT|) . 
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[Figure 5 about here 



In Figure[5l the spectral efficiency of link 1 in the two user MISO IFC setting 
used in Figure[3]and Figure|3]is plotted for perfect and imperfect CSI. The curves 
are found by using the relation Ce{Et/Nog) — Ci{p) for the SNR which solves 
(|40)) . For the plot of Ci{Eij/ Nq) curves, joint MRT beamforming is assumed 
for all values of Eb/No^. Joint MRT is however not Pareto efficient in mid or 
high SNR. From a game theoretic perspective [sO], joint MRT is the dominant 
strategy equilibrium (unique Nash equilibrium) of a strategic game between the 
links [5l|. In other words, if the transmitters are noncooperative, they will 
jointly perform MRT beamforming. Accordingly, the performance plotted in 
Figure [S] is for two noncooperative links. The slopes of the tangents at ^''^^^j^ 

and ^'^.^^j^ correspond to the wideband slopes in Figure 2] with joint MRT for 

perfect and imperfect CSI, respectively. The loss in minimum Eb/No as well as 
wideband slope is apparent due to uncertainty in the channel information. 

5. Conclusions 

We consider a K-nser MISO IFC. The CSI is assumed to be perfect at the 
receivers but imperfect at the transmitters. Channel vector estimates at a trans- 
mitter include channel estimation errors which are assumed to be bounded in 
an elliptical region. The geometry of the uncertainty region associated with 
a channel vector estimate is known at the transmitter. In this setting, robust 
beamforming optimizes worst-case power gains at the receivers. We derive the 
worst-case intended and interference power gains at the receivers and formulate 
accordingly the worst-case achievable rates for the links. Afterwards, we char- 
acterize the robust beamforming vectors necessary to operate at any Pareto 
optimal point in the robust achievable rate region. The efficient beamform- 
ing vectors of each transmitter are found as a solution of a SOCP which can be 
solved efficiently. The spectral efficiency of the multi-link system with imperfect 
channel state information is analyzed in the high and low SNR regime. At high 
SNR, achieving full multiplexing gain with zero forcing transmission using the 
channel vector estimates requires the channel estimation error to reduce linearly 
with the SNR. If the error does not depend on SNR, single-user transmission is 
optimal. In the low SNR regime, it is shown that joint robust maximum ratio 
transmission optimizes the minimum energy per bit for reliable communication. 
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Figure 1: Pareto boundary of a three- user robust rate region at dB SNR. 
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Algorithm 1: Algorithm to calculate the maximum multiplexing gain. 

Input: [Ni,...,Nk] 

1 Sort [Ni , . . . , Nk] in weakly decreasing order to [Ni , . . . , Nk] such that 

Ni > Ni+v, 

2 for k = 1, . . . , K do 

3 ii Nk < k then 

4 m* = k-1; 

5 break; 

Output: m* 



23 



perfect CSI 



-* imperfect CSI with e^.= 1/p 



1/2 



-B — imperfect CSI with e^. = 1/p 

-© — imperfect CSI with e, . = 0.3 

kj 



1/3 



= 2 



S =1 



-20 -10 10 20 30 

signal-to-noise ratio (p) [dB] 



40 



50 



Figure 2: Comparison of the maximum sum rate of two links for different models of channel 
estimation errors e^j, k,j = 1, 2. 
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Figure 3: Plot of minimum energy per bit regions of two links for perfect and imperfect CSI. 
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Figure 5: Plot of the spectral efficiency of link 1 for perfect and imperfect CSI in a two-user 
MISO IFC. 
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